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Abstract 

We investigate the effects of strong magnetic fields upon the gross properties of neutron and 
protoneutron stars. In our calculations, the neutron star matter was approximated by the pure 
neutron matter. Using the lowest order constrained variational approach at zero and finite tem- 
peratures, and employing ^^is potential, we present the effects of strong magnetic fields on the 
gravitational mass, radius, and gravitational redshift of the neutron and protoneutron stars. It is 
found that the equation of state of neutron star becomes stiffer with increase of the magnetic field 
and temperature. This leads to larger values of the maximum mass and radius for the neutron 
stars. 
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I. INTRODUCTION 



Compression of magnetic flux inherited from the progenito r star could form the strong 
magnetic fleld in the interior of a neutron star (Reisenegger 2007l |). Using this point of 
view, Wo ltjer h as predicted a magnetic fleld strength of order 10^^ G for neutron stars 
(Woltjer 1964| ). In the core of high density inhomogeneous gravitatio nally bound neutron 
stars, the magnetic fleld strength can be as large as 10^" G (Ferrer 2010|). In addition, 
considering the formation of a quark core in the high densit y interior of a neu tron star, the 
maximum fleld reaches up to about lO^*' G (Ferrer 2010|; Tatsumi 2000|). According 
to the scalar virial theorem which is based on the Newtonian gravity, the magnetic fleld 
strength is allowed up to 10^^ G in the interior of a magnetar (Lai & Shapiro 199l| ). On 
the other hand, general relativity predicts the allowed maximur n valu e of the neutron star 
magnetic fleld to be about 10^^ — 10^*^ G (Shapiro & Teukolsky 1983| ). By comparing with 
the observational data. Yuan et al. obtained a magnetic fleld strength of order 10^^ G for 
the neutron stars (Yuan & Zhang 1998| ). 

Strong magnetic fleld could have an important influence on the structure of a neutron 
star. Some authors have studied the effects of strong magnetic flelds on the properties of 
neutron stars. Bocquet et al. extended the numerical code for computing the perfect fluid 
rotating neutron stars in general relativity to include the electromagnetic flelds and studie d 
the rapidly rotating neutron stars endowed with magnetic flelds (Bocquet et al. 1995| ). 
The results show that for a magnetic fleld B ~ 10^*^ G, the maximum mass increases by 
13 to 29% (depending upon the equation of state) with respect to the maximum mass 
of non-magnetized stars. Within a relativistic Hartree approach in a simple linear a — 
uj — p model, Chakrabarty et al. studied the gross properties of cold symmetric nuclear 
matter and nuclear ma tter in beta equilibrium under the influence of strong magnetic flelds 

1997|)- They showed that for magnetic flelds B„, = 0, 4.4 x 10^^ 



(Chakrabarty et al. 

and 10^*^ G, the maximum masses are Mmax = B.IOM©, 2.99M0 and 2.91M0, with radii 
^Mmax = 15.02, 14.95, 12.25 km, respectively. Based on two nonlinear a — u models of 
nuclear matter. Yuan et al. considered the p roperties of neutron stars under the influence of 
strong magnetic flelds (Yuan & Zhang 1999| ). They found that the equation of state became 
softer with increase of the magnetic fleld. The results show that for the ZM model, the 



maximum masses are M„ 



I.7OM0 and 1.62Mq for B = 0, 10^" G, with corresponding 



2 



radii R 



M„ 



9.82 and 8.70 km. Furthermore, for the BB model, the maximum masses 



are M^ax = 2.26Mq and 2.O7M0 for B = 0, lO^", with radii Rm^^^ = 12.07, 10.09 km. 
Cardall et al. studied static neutron stars with magnetic fields and a simple class of electr ic 
current distributions consistent with the stationarity requirement (Cardall et al. 200l|). 
It has been demonstrated that the maximum mass of static neutron stars with magnetic 
fields determined by a constant current function is noticeably larger than that attainable 
with uniform rotation and no magnetic field. Within a relativistic field theory, Mao et al. 
considered a neutron-star matter consisting of neutrons, protons and electrons interacting 
through the exchange of a, u and p mesons in the presence of a magn etic field which decreases 
from the center to the surface of a neutron star (Mao et al. [2003| ). It has been found that 
the equation of state becomes stiffer by increasing the magnetic field that led to an increase 
of 40% on the neutron star maximum mass. 

In our previous studies, we have investigated the properties of neu tron stars and pro- 
;oneutr on stars in the ab sence of magnetic field (Bordbar et al. 2006a| ; Bordbar & Hayati 



2006 q : Bordbar et al. [2009 1; Yazdizadeh & Bordbar 2011a |). Recently, we have calcu- 
lated the properties of spi n polar ized neutron matter in the presence of strong m agnet ic fields 
at zero (Bordbar et al. 201 lb| ) and finite temperatures (Bordbar & Rezaei 2012|) using 
LOCV technique employing AVis potential. In the present work, the neutron star matter 
is approximated by the pure neutron matter to investigate the effects of strong magnetic 
fields on the gross properties of neuron stars and protoneutron stars using the equa tions o f 
state of neutron m atter i n the presence of strong magnetic fields (Bordbar et al. 
Bordbar & Rezaei pOlj ). 



2011b) : 



II. NEUTRON STAR STRUCTURE IN THE PRESENCE OF STRONG MAG- 
NETIC FIELDS 



In the present study, we calculate the neutron star and protoneutron star structure using 
the equations of state of cold and hot neutron mat ter iii the presence of strong magnetic 
fields (Bordbar et al. |2011b| : Bo rdbar & Rezaei |2012| ). In our study, we employ AVis 
nuclear potential (Wiringa et al. 1995| ) and use the lowest order constrained variational 



methoc 
et al. 



to ca 



culate the equation of state . For more details, we refer the reader to (Bordbar 



201 lb| : Bordbar & Rezaei |2012| ). Our results for the equation of state of neutron 
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matter in the presence of strong magnetic fields are given in Figs. [UlSl Figs. [U^b) and El^b) 
indicate that for the cold and hot neutron matter, at each density, the pressure increases 
with increase of the magnetic field. This stiffening of the equation of state is due to the 
inclusion of neutron anomalous magnetic moments. In other words, in the presence of high 
magnetic fields, the fraction of polarized neutrons at the equilibrium state increases and 
therefore the degeneracy press ure grows. This is in agreement with the results obtained in 
Refs. (Broderick et al. 2000|; Yue & Shen 2006|). We have found that at low densities, 
the influence of magnetic field on the pressure is negligible. Fig. lU^b) shows that at each 
density, the pressure grows by increasing the temperature. Consequently, for hot neutron 
matter, the equation of state is stiffer compared with the cold one. Fig. El^a) also shows that 
the effect of finite temperature on the equation of state is more significant at high densities. 

The equilibrium configurations could be obtained by solving the general relativistic equa- 
tions of hydrostatic equilibrium, Tolman-Oppenheimer-Volkoff (TOV) (Shapiro & Teukolsky 



[19831) 



dm 
dr 
dP 
dr 



Gm{r)e{r) 



1 + 



F(r) 



A-Kr^Pir) 



1 - 



2Gm(r) 



where e{f) is the energy density, G is the gravitational constant, and 

^Tir''^e(r')dr' 



m(r) 



(2) 



gives the gravitational mass inside a radius r. By selecting a central energy density Ec, under 
the boundary conditions -P(O) = Pc, m(0) = 0, we integrate the TOV equations outwards to 
a radius r = R, at w hich P vanishes. This yields the radius R and mass M = m{R) (Shapiro 
& Teukolsky 1983| ). Gravitational redshift, the criterion for the star compactness, is given 
by 



1, 



(3) 



where R is the radius of the neutron star. In our calculations of neutron star structure, 
for densities greater than 0.05 fm~^, we use the equations of state presented in Figs. [T]- 
ini However, at lower densities, because the magnetic field and finite temperature have 
insignificant e ffects on EoS, we employ the equation of state calculated by Baym et al. 
(Baym et al. 1971| ) for all magnetic fields and temperatures. 



The effects of magnetic fields on the gravitational masses of neutron stars and protoneu- 
tron stars at a temperature about 15 MeV with different central densities are presented in 
Fig. m Obviously, at very low central densities, the gravitational masses are independent of 
the equation of state; but at higher densities, the gravitational mass increases by increasing 
both magnetic field and temperature. The limiting value of neutron star mass (maximum 
mass) also reaches the larger amount when the magnetic field and temperature rise. For a 
cold neutron star at B = 10^^ G, the maximum mass is about 1.17% larger than the cold 
non magnetized one. Considering two stars (a protoneutron star at T = 15 MeV and a cold 
neutron star) in the presence of a magnetic field B = 10^^ G, the protoneutron star maxi- 
mum mass is about 1.16% greater than the cold neutron star. Besides, for a protoneutron 
star at T = 15 MeV and B = 10^^ G, the maximum mass increases about 2.36% compared 
to a cold non magnetized one. These results are due to the stiffening of the equation of state 
(Figs. [H. 

Fig. [5] presents the gravitational mass versus radius (M-R relation) for different magnetic 
fields at zero and finite temperatures. For all magnetic fields and temperatures, the neutron 
star mass decreases by increasing the radius. It is clear from Fig. [5] that for a given radius, 
the gravitational mass increases whenever the equation of state becomes stiffer. We have 
found that the effect of the equation of state on the M-R relation is more significant for the 
neutron stars with smaller radius. 

Fig. [6] shows the gravitational redshift versus the gravitational mass of the neutron star for 
different magnetic fields at zero and finite temperatures. Clearly, the stiffness of the equation 
of state reduces the gravitational redshift. Fig. [6] also indicates that the maximum redshift 
(redshift corresponding to the maximum mass) decreases with the increase of maximum 
mass. For a cold and a hot neutron star at T = 15 MeV with B = 10^^ G, the values of 
maximum redshift are z^"-^ = 0.49, and z"^"-^ = 0.47, respectively. In addition, we have 
found that in the case of a cold neutron star, for magnetic fields B = 0, 5x 10^^, and 10^^ G, 
the values of z^""^ are 0.56, 0.53, and 0.49. Therefore, the maximum surface redshift of our 
calculations, i.e. 34.18% (for a cold neutron star aX B = lO^^G), is lower than the upper 
bound on th e surface redshift for subluminal equation of states, i.e. z^^ = 0.8509 (Haensel 



et al. 



19991). 



Tables 1 and 2 show a summary of our results for the maximum mass and the correspond- 
ing radius predicted for different neutron stars. We have found that the effects of magnetic 
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fields with magnitude B < 10^^ G are almost negligible. Obviously, for cold neutron stars 
as well as protoneutron stars, the maximum mass and the corresponding radius increase by 
increasing the magnetic field. Tables 1 and 2 show that at any magnetic field, the maximum 
mass and the corresponding radius of the protoneuton star are larger than the cold neutron 
star. Therefore, we conclude that the stiffer equation of state leads to a neutron star with 
a larger maximum mass and radius. According to our results, for a cold neutron star, the 
maximum mass can vary between I.GQMq and 1.71Mq, depending on the interior magnetic 
field, but for a protoneutron star with T = 15 MeV, this variation is between I.TOM© and 
I.TSMq. Therefore, the effect of magnetic field on the protoneutron star maximum mass 
is more important than the cold neutron star. Our results for the neutron star maximum 
mass are higher than the observational results from X-ray binaries presented in Table 3. 
Moreover, the study of the statistics of 61 measured masses of neutron stars in binary pulsar 
systems gives a mass average of M = 1.46 ± O.SM© (Zhang et al. [201 1| ). Their results 
indicate that the mass average of the more rapidly rotating millisecond pulsars (MSPs) is 
M = 1.57±O.35M0. In the present work the values of the protoneutron star radius at highe r 



magnetic fields are near the values obtained using MR relationships (Zhang et al. [200 7l |) 
which shows the neutron star radius varies in the range of 10 — 20 km. We have also found 
that the effect of magnetic field on the radius of the protoneutron star is less important than 
the cold neutron star. 



III. SUMMARY AND CONCLUSION 



Different properties of the neutron star and protoneutron star structure have been investi- 
gated using the equation of state of neutron matter in the presence of strong magnetic fields. 
In our calculations, we have employed the lowest order constrained variational method and 
applied AVis potential to find the equation of state at zero and finite temperature in the 
presence of strong magnetic fields. Our results show that the stiffer equation of state at 
higher magnetic fields and larger values of temperatures lead to the higher values for the 
maximum mass and radius. For the maximum value of the magnetic field considered in this 
study, i.e. 10^^ G, the maximum masses of a cold neutron star and protoneutron star at 
T = 15 MeV are I.TIM© and I.TSMq, respectively. The corresponding radii are also 9.16 
and 9.22 km. Our results indicate that the effects of magnetic field on the maximum mass 
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of the protoneutron stars are more important than cold neutron stars, while the effects of 
magnetic fields are more visible on the radius of cold neutron stars. It has been shown that 
the effects of the equation of state on the M-R relation are more important for neutron 
stars with smaller radii. Our calculations also demonstrate that the maximum value of the 
gravitational surface redshift decreases by increasing the neutron star maximum mass. 
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TABLE I: Maximum gravitational mass, Mmax, and the corresponding radius, RMmaxi obtained 
for different values of magnetic field, i?, at T = MeV . 



B{G) 


M^ax{MQ) 







1.69 


8.59 


5 X 101*^ 


1.70 


8.73 


1019 


1.71 


9.16 



TABLE 11: Same as Table 1 but at T = IhMeV. 



B{G) 


M^ax{MQ) 







1.70 


8.70 


5 X 10^^ 


1.71 


8.83 


1019 


1.73 


9.22 



TABLE TIL Measured masses of neutron stars in X-ray binaries. 



System 


M{Mq) 


References 


SMC X-1 
Cen X-3 
LMC X-4 
V395 CAR/2S 0921C630 


1.05 ±0.09 
1.24 ±0.24 
1.31 ±0.14 
1.44 ±0.10 


(van der Meer et al. [2005]) 
(van der Meer et al. [2005]) 
(van der Meer et al. [2005]) 
(Steeghs et al. [2007]) 
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FIG. 1: (a) Pressure, P, versus energy density, e, for the cases B = G (solid curve), B = 5x IQ-*^^ G 
(dashed curve) and B = 10^^ G (dashdot curve) at a fixed value of the temperature, T = MeV. 
(b) Same as in the top panel but for a different range of density. 



10 



J \ \ I \ \ \ \ I \ \ \ \ I \ \ \ \ I \ \ \ ^1 \ \ \ \ I \ 1 

5 10 15 20 25 30 

8 (10^' g/cm^) 




e(10'' g/cm^) 



FIG. 2: (a) Pressure, P, versus energy density, e, for the cases B = G (solid curve), B = 5x IQ-*^^ G 
(dashed curve) and B = 10^^ q (dashdot curve) at a fixed value of the temperature, T = 15 MeV. 
(b) Same as in the top panel but for a diff^erent range of density. 
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FIG. 3: (a) Pressure, P, versus energy density, e, for the cases T = MeV (solid curve) and 
T = 15 MeV (dashed curve) at a fixed value of the magnetic field, B = 5 x 10^^ G. (b) Same as 
in the top panel but for a different range of density. 
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FIG. 4: (a) Gravitational mass of neutron star (in units of the solar mass, Mq) versus central 
energy density, Ec, at T = MeV. All curves correspond to those of Fig. [TJ (b) Same as (a) but 
at T = 15 MeV. All curves correspond to those of Fig. [21 (c) Gravitational mass of neutron star 
(in units of the solar mass, Mq) versus central energy density, Ec, at B = 5 x 10^® G. All curves 
correspond to those of Fig. [3l 
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FIG. 5: (a) Mass-radius relation at T = MeV. All curves correspond to those of Fig. [TJ (b) 
Same as (a) but at T = 15 MeV . All curves correspond to those of Fig. [2j (c) Mass-radius relation 
at i? = 5 X 10^*^ G. All curves correspond to those of Fig. [3l 
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FIG. 6: Gravitational redshift, Z^, vs. total mass for neutron stars at T = MeV. All curves 
correspond to those of Fig. [TJ (b) Same as (a) but at T = 15 MeV. All curves correspond to those 
of Fig. [2j (c) Gravitational redshift, Zs, vs. total mass for neutron stars at = 5 x 10^^ G. All 
curves correspond to those of Fig. [3l 
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